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Abstract 

We evaluate differential distributions for exclusive scalar /o(1500) production for pp — ► iVi^/o 
(FAIR@GSI) and pp — > ppfo (J-PARC@Tokai). Both QCD diffractive and pion-pion meson ex- 
change current (MEC) components are included. Rather large cross sections are predicted. The 
pion-pion component, never discussed in the literature, dominates close to the threshold while the 
diffractive component takes over for larger energies. The diffractive component is calculated based 
on two-gluon impact factors as well as in the framework of Khoze-Martin-Ryskin approach pro- 
posed for diffractive Higgs boson production. Different unintegrated gluon distribution functions 
(UGDFs) from the literature are used. The production of /o(1500) close to threshold could limit 
the so-called irNN form factor in the region of larger pion virtualities. 
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I. INTRODUCTION 



Many theoretical calculations, including lattice QCD, predicted existence of glueballs 
(particles dominantly made of gluons) with masses M > 1.5 GeV. No one of them was up to 
now unambigously identified. The lowest mass meson considered as a glueball candidate is a 
scalar / (1500) [5J discovered by the Crystall Barrel Collaboration in proton-antiproton an- 
nihilation [l|. The branching fractions are consistent with the dominant glueball component 
0]. It was next observed by the WA102 collaboration in central production in proton- proton 
collisions in two-pion 0] and four-pion [ij decay channels at -Js ~ 30 GeV 1 . Close and Kirk 
proposed a phenomenological model of central exclusive /o(1500) production. In their 
language the pomerons (transverse and longitudinal) are the effective (phenomenological) 
degrees of freedom [7|. The Close-Kirk amplitude was parametrized as 

M(t 1 ,t 2 ,<f)')=aTexp(^(t 1 +t2)\ + a L ^^-exp\^-(t 1 +t 2 )\ cos(0') • (1-1) 

In their approach there is no explicit /o(1500)-rapidity dependence of the corresponding am- 
plitude. Since the parameters were rather fitted to the not-normalized WA102 data 0] no ab- 
solute normalization can be obtained within this approach. Furthermore the parametrization 
is not giving energy dependence of the cross section, so predictions for other (not-measured) 
energies are not possible. In the present paper we will concentrate rather on a QCD-inspired 
approach. It provides absolute normalization, energy dependence and dependence on meson 
rapidity (or equivalently on xp of the meson). 

The nature of the /o(1500) meson still remains rather unclear. New large-scale devices 
being completed (J-PARC at Tokai) or planned in the future (FAIR at GSI) open a new 
possibility to study the production of /o(1500) in more details. 

In the present analysis we shall concentrate on exclusive production of scalar /o(1500) in 
the following reactions: 

p + p^p + /o(1500)+p, 

p+p->p + /o(1500)+p, (1.2) 
p + p -> n + / (1500) + n , 

While the first process can be measured at J-PARC, the latter two reactions could be mea- 
sured by the PANDA Collaboration at the new complex FAIR planned in GSI Darmstadt. 
The combination of these processes could shed more light on the mechanism of /o(1500) 
production as well as on its nature. 

If /o(1500) is a glueball (or has a strong glueball component (9)]) then the mechanism 
shown in Fig. [T]may be important, at least in the high-energy regime. This mechanism is 
often considered as the dominant mechanism of exclusive Higgs boson ^ and Xc(0 + ) meson 
[ll| production at high energies. There is a hope to measure these processes at LHC in 
some future when forward detectors will be completed. At intermediate energies the same 
mechanism is, however, not able to explain large cross section for exclusive 7/ production 
l~ol | as measured by the WA102 collaboration. Explanation of this fact is not clear to us in 



the moment. 



No absolute normalization of the corresponding experimental cross section was available. Only two-pion 
or four-pion invariant mass spectra were discussed. 
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FIG. 1: The sketch of the bare QCD mechanism for diffractive production of the glueball. The 
kinematical variables are shown in addition. 



At lower energies (y/s < 20 GeV) other processes may become important as well. Since 
the two-pion channel is one of the dominant decay channels of /q(1500) (34.9 ± 2.3 %) 
22| one may expect the two-pion fusion (see Figf5]to be one of the dominant mechanisms of 



exclusive / (1500) production at the FAIR energies. The two-pion fusion can be also relative 
reliably calculated in the framework of meson exchange theory. The pion coupling to the 
nucleon is well known [l5|. The ttNN form factor for larger pion virtualities is somewhat less 
known. This may limit our predictions close to the threshold, where rather large virtualities 
are involved due to specific kinematics. At largest HESR (antiproton ring) energy, as will be 
discussed in the present paper, this is no longer a limiting factor as average pion virtualities 
are rather small. 
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FIG. 2: The sketch of the pion-pion MEC mechanism. Form factors appearing in different vertices 
and kinematical variables are shown explicitly. 
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II. EXCLUSIVE PROCESSES 



A. Cross section and phase space 

The cross section for a general 3-body reaction pp — > pp / (1500) can be written as 

da pp ^ ppM = —^=2===]M^- d*PS . (2.1) 

Above m is the mass of the nucleon. 

The three-body phase space volume element reads 

d 3 PS = y d 3p M . (2vr )^4 ( + / _ / _ P s (2.2) 

2^(27r) 3 2^(2vr) 3 2 J B M (27r) 3 V ; yF y Fl F2 ' V ; 

At high energies and small momentum transfers the phase space volume element can be 



written as 14 



d3pS ~ ^dtxdhdttdbdcj) 6 - ^i)(l - 6) - M 2 ) , (2.3) 

where £i, £2 are longitudinal momentum fractions carried by outgoing protons with respect to 
their parent protons and the relative angle between outgoing protons $ e (0, 2tt). Changing 
variables (^1,^2) - ► (%f,M 2 ) one gets 

d 3 PS « -^-rdhdh — dXF d$ . (2.4) 

2 8 ^ 4 V4 + 4(M2 + |P M ,|2)/ S ^ ' 

The high-energy formulas ( 12.31) and (12.41) break close to the meson production thresh- 
old. Then exact phase space formula ( 12.21) must be taken and another choice of variables 
is more appropriate. We choose transverse momenta of the outgoing nucleons (Pi t ,P2t)> az ~ 
imuthal angle between outgoing nucleons (0) and rapidity of the meson (y) as independent 
kinematically complete variables. Then the cross section can be calculated as: 



7 rr-X( , \ I \M(p lti P2U^y)\ 2 27T 1 1 1 A A AAA (O K\ 

da=yj {pit,P2tA,y)\k ~ -1=, == . > 5 , OT? , -pitP2tdpitdp 2 td(j)dy , (2.5) 

^ 2^s(s — Am 1 ) lh, x 2ti 2 I 

where k denotes symbolically discrete solutions of the set of equations for p' lz and p' 2z : 

Vs-Em = ^m\ t +p' 2 z + y/ml t +pg , , 2 ^ 

-PMz = p' lz + p' 2z , 

where m\ t and rri2t are transverse masses of outgoing nucleons. The solutions of 
Eq. (12.61) depend on the values of integration variables: p' lz = p' lz (p' lt ,p' 2t ,4>,y) and p' 2z = 
p'lJ^P'xtiPui 4>i v)- The extra jacobian reads: 



Jk 



P'lM P f 2z(k) 



V m it + p[z( k ) 2 V m it+P2z( k y 



(2.7) 



In the limit of high energies and central production, i.e. p' lz ^> (very forward nucleonl), 
—p' 2z ^> (very backward nucleon2) the jacobian becomes a constant J — > ^. 

The matrix element depends on the process and is a function of kinematical variables. 
The mechanism of the exclusive production of /o(1500) close to the threshold is not known. 
We shall address this issue here. Therefore different mechanisms will be considered and the 
corresponding cross sections will be calculated. 
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B. Diffractive QCD amplitude 



According to Khoze-Martin-Ryskin approach (KMR) |8j], we write the amplitude of ex- 
clusive double diffractive colour singlet production pp — > pp / (1500) as 



M 



a 9' 



s 2 1 5, 
2' 71 2iV~ 2 



C\C2 rx. 



fg,l X 'li <?0,t> 1l,t' ^)fgji( X ^ X 2' 3o,t' ?2,t' 



off. 



d 2 qo,tV C j 



qlt Qlt qlt 



(2.8) 



The normalization of this amplitude differs from the KMR one [8J by the factor s/ 2 and 
coincides with the normalization in our previous work on exclusive ^'-production [10] . The 
amplitude is averaged over the colour indices and over two transverse polarisations of the 
incoming gluons The bare amplitude above is subjected to absorption corrections which 
depend on collision energy (the bigger the energy, the bigger the absorption corrections). 
We shall discuss this issue shortly when presenting our results. 

The vertex factor Vj 1C2 = Vj lC2 (qf t , q 2t , P%f t ) in expression (12.81) describes the coupling of 
two virtual gluons to /o(1500) meson. Recently the vertex was obtained for off-shell values 
of qxj and q 2jt in the case of x c (0) exclusive production [llj]. An almost alternative way 
to describe the vertex is to express it via partial decay width T(M — > gg). 2 The latter 
(approximate) method can be used also for glueball production. 

In the original Khoze-Martin-Ryskin (KMR) approach [8j the amplitude is written as 



M = N 



rfV t P[/o(150 0)] 

q\ t q\A 



f« MR (x u x[, Ql t , ^ h)ff MR { X2 , x' 2 , Qi t , ^ t 2 



(2.9) 



where only one transverse momentum is taken into account somewhat arbitrarily as 



Qlt = mm i9o,t> 9i,J > Ql,t = mm i9o,i> qlt) , 

and the normalization factor iV can be written in terms of the /o(1500) - 
(see below). 

In the KMR approach the large meson mass approximation M ^> |qi^| 
so the gluon virtualities are neglected in the vertex factor 



P[/o(1500)] ~ ( 9l>t g 2)t ) = (qo,t+p[, t )(qo,t 
The KMR UGDFs are written in the factorized form: 



P2, 



f« MR (x, x', Ql /i 2 ; t) = f? MR (x, x 1 , Ql /x 2 ) exp(M) 



(2.10) 
gg decay width 

\q,2,t\ is adopted, 
(2.11) 

(2.12) 



with bo = 2 GeV~ 2 [8] . In our approach we use somewhat different parametrization of the 
t-dependent isoscalar form factors. 

Please note that the KMR and our (general) skewed UGDFs have different number of 
arguments. In the KMR approach there is only one effective gluon transverse momentum (see 
Eq. (l2.10p ) compared to two idependent transverse momenta in general case (see Eq. fl2.16l) ). 

The KMR skewed distributions are given in terms of conventional integrated densities g 
and the so-called Sudakov form factor T as follows: 



f« MR (x,x>,Ql^) 







d\nQ 2 t 



T(Qlfi 2 )xg(x,Ql 



(2.13) 



2 The last value is not so well known. We shall take T(M — > gg) = T l £} . This will give us un upper estimate. 
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The square root here was taken using arguments that only survival probability for hard 
gluons is relevant. It is not so-obvious if this approximation is reliable for light meson 
production. The factor R g in the KMR approach approximately accounts for the single 
logQ 2 skewed effect (sf. Please note also that in contrast to our approach the skewed 
KMR UGDF does not explicitly depend on x' (assuming x' <C x <C 1). Usually this 
factor is estimated to be 1.3—1.5. In our evaluations here we take it to be equal 1 to avoid 
further uncertainties. Following now the KMR notations we write the total amplitude ( 12. 8|) 
(averaged over colour and polarisation states of incoming gluons) in the limit M ^> qi )t , g 2 ,t 
as 

M = Atx 2 - / d 2 q 0it P [f 1500 : ' 2 2 ; ! , (2.14) 

where the normalization constant is 



a _ 647rr(/ (1500) -> gg) 
A ~ (N* - 1)M3 • (2 - 15) 

In addition to the standard KMR approach we could use other off-diagonal distributions 
(for details and a discussion see [Io|, [Til] ) - In the present work we shall use a few sets of 
unintegrated gluon distributions which aim at the description of phenomena where small 
gluon transverse momenta are involved. Some details concerning the distributions can be 
found in Ref. T3l- We shall follow the notation there. 



In the general case we do not know off-diagonal UGDFs very well. In [10|, [ll| we have 
proposed a prescription how to calculate the off-diagonal UGDFs: 

fl" = \lf?\^qlt^l)-ti\^ql t ^ 2 ) ■ F^h) , (2.16) 

where and Fi(t 2 ) are isoscalar nucleon form factors. They can be parametrized as 

(0) 

Am 2 p ~ 2 - 79t i,2 

-^1(^1,2) = j- A — 2 7 — \7i ' / n71 \ 2 ■ ( 2 -^7) 

(Am 2 p - ti, 2 )(l - ti, 2 /071) 2 

Above t\ and t 2 are total four-momentum transfers in the first and second proton line, 
respectively. While in the emission line the choice of the scale is rather natural, there is no 
so-clear situation for the second screening-gluon exchange flpj ] . 

Even at intermediate energies (W = 10-50 GeV) typical relatively small 

(~ 0.01). However, characteristic x\,x 2 ~ M/ /-y/s are not too small (typically > 10 _1 ). 
Therefore here we cannot use the small-x models of UGDFs. In the latter case a Gaussian 
smearing of the collinear distribution seems a reasonable solution: 

^ auss (x, kl /4) = xg coll (x, /4) • f GauS s(k 2 ; a ) , (2.18) 

where g coU (x , fi 2 F ) are standard collinear (integrated) gluon distribution and fcaussik 2 ; a ) is 
a Gaussian two-dimensional function 

f G auss(k 2 , <7 ) = - — 2- exp (-k 2 /2a 2 ) /it . (2.19) 

Z7TCT q 
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Above o"o is a free parameter which one can expect to be of the order of 1 GeV. Based on 
our experience in 10J we expect strong sensitivity to the actual value of the parameter &q. 
Summarizing, a following prescription for the off-diagonal UGDF seems reasonable: 



f(x jX ',klk t \t) = Jf sm all- X (x',K 2 )fGauss(x,kl^) ■ F(t) , (2.20) 



where f S maU-x{x' , k t 2 ) is one of the typical small- x UGDFs (see e.g.(l2j). So exemplary 
combinations are: KL ® Gauss, BFKL cg> Gauss, GBW <g> Gauss (for notation see |12|). 
The natural choice of the scale is \j? = M? . This relatively low scale is possible with the 



GRV-type of PDF parametrization [13J. We shall call (12.201) a "mixed prescription" for 
brevity. 



C. Two-gluon impact factor approach for subasymptotic energies 

The amplitude in the previous section, written in terms of off-diagonal UGDFs, was 
constructed for large energies. The smaller the energy the shorter the QCD ladder. It is 
not obvious how to extrapolate the diffractive amplitude down to lower (close-to-threshold) 
energies. Here we present slightly different method which seems more adequate at lower 
energies. 

At not too large energies the amplitude of elastic scattering can be written as amplitude 



for two-gluon exchange [2Cj, l21[ 

awm) = / dX a - (k " )a - {k - ) 3J m^^0 1 ■ (2 - 21) 

In analogy to dipole-dipole or pion-pion scattering (see e.g. [21] ) the impact factor can be 
parametrized as: 

A 2 A 2 
F(k "' ^ = A 2 + (k u + k M )* - A 2 + (k u - k 2t)2 • (Z22) 

At high energy the net four-momentum transfer: t = — (k u + k 2 t) 2 . A in Eq. (12.211) is a 
free parameter which can be adjusted to elastic scattering. For our rough estimate we take 
A = m p . 

Generalizing, the amplitude for exclusive /o(1500) production can be written as the am- 
plitude for three-gluon exchange shown in Figj3j 



M pp - m> f (i5oo)(s,y,ti,h,<f) = is " C N 2 J rf2fc °* 

(a s {kl t )a s (kl t )) 1/2 {a s (kl t )a s (k 2 2t )) 
3F(k ot , k u )3F(k ot , k 2t 



(kl + tf)(kl t + tf)(kl + ti 



N 2 

1/2 M4H(4)) 1/2 

^gs^/o (1500) (kit, k 2t ) • (2.23) 



At high energy and y « the four-momentum transfers can be calculated as: 
t x = -(ko* + k 14 ) 2 , h = -(kw - k 2t ) 2 . 
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FIG. 3: The sketch of the two-gluon impact factor approach. Some kinematical variables are shown 
explicitly. 

At low energy and/or y ^ the kinematics is slightly more complicated. Let us define 
effective four-vector transfers: 

91 = (p'i -Pi) = (qio,qix,qi y ,qiz) , 

92 = (p' 2 -P2) = (920,922, 42?,, fe) • (2.24) 

Then t x = q\ = q\ x + q\ t and t 2 = q\ = q 2l + q 2t . Close to threshold the longitudinal 
components qf t = qf — q\ z <C and q 2l = q 20 — q\ z -C 0. Then the amplitude ( 12. 23ft must 
be corrected. Then also four-vectors of exchanged gluons (k , k x and k 2 ) cannot be purely 
transverse and longitudinal components must be included as well. To estimate the effect 
we use formula ( 12.231) 3 but modify the transferred four momenta of gluons entering the 
9*9* ~^ /o(1500) production vertex: 

fci = (0, kit, 0) ->• (gio, k lt , q lz ) , 

k 2 = (0,k 2t ,0) -> (g 2 o,k 2i ,g 22 ) (2.25) 

and leave ko purely transverse. This procedure is a bit arbitrary but comparing results 
obtained with formula (I2.23j) with that from the formula with modified four-momenta would 
allow to estimate related uncertainties. 

We write the vertex function gg — > / (1500) in the following tensorial form 4 : 

V(h, k 2 ,p M ) = C Ml500 )->gg g^k v 2 . (2.26) 

The normalization factor is obtained from the decay of /o(1500) into two soft gluons: 

|C/o(1500)^3 3 | 2 = M 3 l N 2 _ jT r /o(1500)-* flfl • (2.27) 



3 It would be more appropriate to calculate in this case a four-dimensional integral instead of the two- 
dimensional one. 

4 In general, another tensorial forms are also possible. This may depend on the structure of the considered 
meson. 



S 



Of course the partial decay width is limited from above: 



r/ (1500)^gg < ^tot ■ (2.28) 

The amplitudes discussed here involve transverse momenta in the infra-red region. Then 
a prescription how to extend the perturbative a s (fc|) dependence to a nonperturbative region 
of small gluon virtualities is unavoidable. In the following a s (kf) is obtained from an analytic 
freezing proposed by Shirkov and Solovtsev Ijl • 



D. Pion-pion MEC amplitude 

It is straightforward to evaluate the pion-pion meson exchange current contribution shown 
in Figj2j If we assume the ij*, type coupling of the pion to the nucleon then the Born 
amplitude reads: 

\M?=\ m + m) (K + m) (^_ + _ + 2 P' )( g + m) )] ■ 2 

77 ZIry2 p -KNN\ t i) ■ 1^/0(1500)^1 V-K^foiisw^M) ■ 77 —^^■kNnK^V 

[(E 2 + m)(E>+m) ( —if- + - 2pa ; p )]-2 

\(E 2 + m) 2 (E 2 + m) 2 (E 2 + m)(E 2 + m) J 

(2.29) 

In the formula above m is the mass of the nucleon, E±, E 2 and E[, E' 2 are energies of initial 
and outgoing nucleons, Pi,P2 and p'^p^ are corresponding three- momenta and is the 
pion mass. The factor q^nn is the familiar pion nucleon coupling constant which is precisely 

2 

known {^f- = 13.5 - 14.6). The isospin factor equals 1 for the 7r°7r° fusion and equals 2 
for the 7r + 7r~ fusion. In the case of proton-proton collisions only the 7r°7r° fusion is allowed 
while in the case of proton-antiproton collisions both tc°tt and 7r + 7r~ MEC are possible. 
In the case of central heavy meson production rather large transverse momenta squared 
ti and t 2 are involved and one has to include extended nature of the particles involved in 
corresponding vertices. This is incorporated via -FViViV^i) or F 1TNN {t 2 ) vertex form factors. 
The influence of the t-dependence of the form factors will be discussed in the result section. 
In the meson exchange approach [l6[ they are parametrized in the monopole form as 

A 2 _ ^2 

F*Mt) = ~J^~f ■ ( 2 - 3 °) 



A typical values are A = 1.2—1.4 GeV [16J. The Gottfried Sum Rule violation prefers smaller 



A « 0.8 GeV 17 



The normalization constant \C\ 2 in (12.291) can be calculated from the partial decay width 

as 

\ r |2 _ 8tt 2M| o r /o (iboo)^^ 
1^/0(1500)^7™! — 1 = > 

\j M l ~ 4 < 

where rj ( 1500 )_ >7r o 7r o = 0.109 • BR(f (1500) — > nn) ■ 0.5 GeV. The branching ratio is 
B R(f o(1500) — > 7T7r) = 0.349 (22|. The off-shellness of pions is also included for the 
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7T7T — > / (1500) transition through the extra K r7r ^/ (i 50 o)(ii, ^2) form factor which we take in 
the factorized form: 



AL f -m! A 



^TTTr-^/o (1500) ^1^2; — -^2 ~l A 1 ' 

It is normalized to unity when both pions are on mass shell 

V(t 1 =ml,t 2 = ml) = l. (2.33) 
In the present calculation we shall take A^^, = 1.0 GeV. 



III. RESULTS 

A. Gluonic QCD mechanisms 

Let us start with the QCD mechanism relevant at higher energies. We wish to present 
differential distributions in xp, ti or i 2 and relative azimuthal angle 0. In the following 
we shall assume: r / (i5oo)->gg = ^f^oo) • This assumption means that our differential 
distributions mean upper value of the cross section. If the fractional branching ratio is 
known, our results should be multiplied by its value. 

In FigJH we show as example distribution in Feynman xp for Kharzeev- Levin UGDF 
(solid) and the mixed distribution KL ® Gaussian (dashed) for several values of collision 
energy in the interval W = 10 - 50 GeV. In general, the higher collision energy the larger 
cross section. With the rise of the initial energy the cross section becomes peaked more and 
more at Xp ~ 0. The mixed UGDF produces slightly broader distribution in xp. 

In Fig J5] we present corresponding distributions in t — t\ = t 2 . The slope depends on 
UGDF used, but for a given UGDF is almost energy independent. 

Finally we present corresponding distributions in relative azimuthal angle between outgo- 
ing protons or proton and antiproton 5 . These distributions have maximum when outgoing 
nucleons are back-to-back. Again the shape seems to be only weekly energy dependent. 



B. Gluonic versus pion-pion mechanism 

What about the pion-pion fusion mechanism? Can it dominate over the gluonic mecha- 
nism discussed in the previous subsection? In FigJ7]we show the integrated cross section for 
the exclusive / (1500) elastic production 

pp^p/ (1500)p (3.1) 

and for double charge exchange reaction 

VV~^ n/o(1500)n • (3.2) 

The thick solid line represents the pion-pion component calculated with monopole vertex 
form factors (12.301) with A = 0.8 GeV (lower) and A = 1.2 GeV (upper). The difference 



The QCD gluonic mechanism is of course charge independent. 
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FIG. 4: The distribution of / (1500) in Feynman x F for W = 10, 20, 30, 40, 50 GeV. In this 
calculation the Kharzeev-Levin UGDF (solid line) and the mixed distribution KL® Gauss (dashed 
line) were used. 

between the lower and upper curves represents uncertainties on the pion-pion compenent. 
The pion-pion contribution grows quickly from the threshold, takes maximum at W ~ 6-7 
GeV and then slowly drops with increasing energy. The gluonic contribution calculated 
with unintegrated gluon distributions drops with decreasing energy towards the kinematical 
threshold and seems to be about order of magnitude smaller than the pion-pion component 
at W = 10 GeV. We show the result with Kharzeev-Levin UGDF (dashed line) which 
includes gluon saturation effects relevant for small-x, Kimber-Martin-Ryskin UGDF (dotted 
line) used for the exclusive production of the Higgs boson and the result with the "mixed 
prescription" (KL ® Gaussian) for different values of the cr parameter: 0.5 GeV (upper thin 
solid line), 1.0 GeV (lower thin solid line). In the latter case results rather strongly depend 
on the value of the smearing parameter. 

We calculate the gluonic contribution down to W = 10 GeV. Extrapolating the gluonic 
component to even lower energies in terms of UGDFs seems rather unsure. At lower energies 
the two-gluon impact factor approach seems more relevant. The impact factor approach 
result is even order of magnitude smaller than that calculated in the KMR approach (see 
lowest dash-dotted (red on-line) line in Fig. [7]), so it seems that the diffractive contribution 
is rather negligible at the FAIR energies. 

Our calculation suggests that quite different energy dependence of the cross section may 
be expected in elastic and charge-exchange channels. Experimental studies at FAIR and 
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FIG. 5: Distribution in t = t x = t 2 for Kharzeev-Levin UGDF for W = 10, 20, 30, 40, 50 GeV. 
The notation here is the same as in FigUl 

J-PARC could shed more light on the glueball production mechanism. 
C. Predictions for PANDA at HESR 

Let us concentrate now on pp collisions at energies relevant for future experiments at 
HESR at the FAIR facilite in GSI. Here the pion-pion MEC (see Figj2]) seems to be the 
dominant mechanism, especially for the charge exchange reaction pp —>■ nnf (1500). 

In FigJH]we show average values of t\ (or £2) for the two-pion MEC as a function of the 
center of mass energy. Close to threshold W = 2mAr+m/ (i5oo) the transferred four-momenta 
squared are the biggest, of the order of about 1.5 GeV 2 . The bigger energy the smaller the 
transferred four-momenta squared. Therefore experiments close to threshold open a unique 
possibility to study physics of large transferred four-momenta squared at relatively small 
energies. This is a quite new region, which was not studied so far in the literature. 

The maximal energy planned for HESR is ^/s = 5.5 GeV. At this energy the phase space 
is still very limited. In Figj9]we show rapidity distribution of /o(1500). For comparison the 
rapidity of incoming antiproton and proton is 1.72 and -1.72, respectively. This means that 
in the center-of-mass system the glueball is produced at midrapidities, on average between 
rapidities of outgoing nucleons. 

In FigHU] we show transverse momentum distribution of neutrons or antineutrons pro- 
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FIG. 6: Distribution in relative azimuthal angle for different UGDFs for W = 10, 20, 30, 40, 50 
GeV. The notation here is the same as in FigUJ 

duced in the reaction pp — > nnf (1500). The distribution depends on the tcNN form factors 
F n NN(h) an d F nNN (t 2 ) in formula (j2.29j) . 

In Fig[IT] we show azimuthal angle correlation between outgoing hadrons (in this case 
neutron and antineutron) . The preference for back-to-back configurations is caused merely 
by the limitations of the phase space close to the threshold. This correlation vanishes in the 
limit of infinite energy. In practice far from the threshold the distribution becomes almost 
constant in azimuth. This has to be contrasted with similar distributions for pomeron- 
pomeron fusion shown in FigJH] which are clearly peaked for the back-to-back configurations. 
Therefore a deviation from the constant distribution in relative azimuthal angle for the 
highest HESR energy of W = 5.5 GeV for pp — > pf (1500)p can be a signal of the gluon 
induced processes. It is not well understood what happens with the gluon induced diffractive 
processes when going down to intermediate (W = 5-10 GeV) energies. A future experiment 
performed by the PANDA collaboration could bring new insights into this issue. This would 
be also a signal that the /o(1500) state has a considerable glueball component. 

Up to now we have neglected interference between pion-pion and pomeron-pomeron con- 
tributions (for the same final channel). This effect may be potentially important when both 
components are of the same order of magnitude. While the pomeron-pomeron contribution 
is dominantly nucleon helicity preserving the situation for pion-pion fusion is more compli- 
cated. In the latter case we define 4 classes of contributions with respect to the nucleon 
helicities: cc (both helicity conserved), cf (first conserved, second flipped), fc (first flipped, 
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FIG. 7: The integrated cross section as a function of the center of mass energy for pp — ► pp/o(1500) 
(left panel) and pp — ► nn/o(1500) (right panel) reactions. The thick solid lines are for pion-pion 
MEC contribution (A = 0.8, 1.2 GeV), the dashed line is for QCD diffractive contribution obtained 
with the Kharzeev-Levin UGDF, the dotted line for the KMR approach and the thin solid lines 
(blue on-line) are for "mixed" UGDF (KL (g) Gaussian) with o"o = 0.5, 1 GeV. The dash-dotted 
line represents the two-gluon impact factor result. 

second conserved) and // (both helicities nipped). The corresponding ratios of individual 
contributions to the sum of all contributions are shown in FigJTiZl In practice, only the cc 
mi contribution may potentially interfere with the gluonic one. From the figure one can 
conclude that this can happen only when both transverse momenta of the final nucleons are 
small. We shall leave numerical studies of the interference effect for future investigations, 
when experimental details of such measurements will be better known; but already now one 
can expect them to be rather small. 



IV. DISCUSSION AND CONCLUSIONS 

We have estimated the cross section for exclusive /o(1500) production not far from the 
threshold. We have included both gluon induced diffractive mechanism and the pion-pion ex- 
change contributions. The first was obtained by extrapolating down the cross section in the 
Khoze-Martin-Ryskin approach with unintegrated gluon distributions from the literature as 
well as using two-gluon impact factor approach. A rather large uncertainties are associated 
with the diffractive component. The calculation of MEC contribution requires introducing 
extra vertex form factors which are not extremely well constraint. This is especially impor- 
tant close to the threshold where rather large \ti\ and 1^1 are involved. The cross section 
for energies close to the threshold is very sensitive to the functional form and parameters of 
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FIG. 8: Average value of < t% >=< t<i > as a function of the center-of-mass collision energy for 
the two-pion exchange mechanism. 

vertex form factor. Therefore a measurement of /o(1500) close to its production threshold 
could limit the so-called nNN form factors in the region of exchanged four-momenta never 
tested before. 

We predict the dominance of the pion-pion contribution close to the threshold and diffrac- 
tive component far from the threshold. Taking into account rather large uncertainties these 
predictions should be taken with some grain of salt. Clearly an experimental program is 
required to disantagle the reaction mechanism. 

Disantangling the mechanism of the exclusive /o(1500) production not far from the me- 
son production threshold would require study of the pp — > pp/o(1500), pp — > nnf Q (1500) 
processes with PANDA detector at FAIR and pp — > ppfo(l500) reaction at J-PARC. In the 
case the gluonic mechanisms are small and the pion exchange mechanism is a dominant 
process one expects: a{pp — > nn/o(1500)) = 4 x a{pp — > pp/o(1500)). On the other hand 
if the gluonic components dominate over MEC components a(pp — > ppj o(1500)) > a(pp — > 
nn/ (1500)). 

Therefore a careful studies of different final channels at FAIR and J-PARC could help to 
shed light on coupling of (nonperturbative) gluons to / (1500) and therefore would give a 
new hint on its nature. Such studies are not easy at all as in the im decay channel one expects 
a large continuum. This continuum requires a better theoretical estimate. A partial wave tttc 
analysis may be unavoidable in this context. The two-pion continuum will be studied in our 
future work. A smaller continuum may be expected in the KK or four-pion /o(1500) decay 
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FIG. 9: rapidity distribution of /o(1500) produced in the reaction pp — » nn/o(1500) for W = 3.5, 
4.0, 4.5, 5.0, 5.5 GeV (maximal HESR energy). In this calculation A = 0.8 GeV. 

channel. This requires, however, a good geometrical (full solid angle) coverage and high 
registration efficiencies. PANDA detector seems to fullfil these requirements, but planning 
real experiment requires a dedicated Monte Carlo simulation of the apparatus. 
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FIG. 12: Helicity decomposition of the cross section on the (pit,P2t) plane for W = 10 GeV. R cc 
(upper left), Rff (upper right), R c f (lower left), Rf c (lower right). The standard nucleon dipole 
form factor was used in this calculation. 



19 



